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as Tycho 651672 and GSC632 1-0 1030) was successfully observed 
with the Hubble Space Telescope (HST) in daylight over the 
northwest Pacific Ocean; Astrometer 3 of the Fine Guidance Sensors 
(FGS) was used to record the event 5 . Details of the HST data are 
given in Table 1 along with information about other observing 
stations; the immersion and emersion data (disappearance and 
reappearance of the star) are shown in Fig. 1. A central flash (the 
focusing of light rays by Triton’s atmosphere 6,7 ) was recorded, but 
will be presented and analysed elsewhere. 

We modelled the HST light curve with a standard small -planet 
model that allows for a power-law temperature gradient 8 (Table 2). 
The background from dark counts and Triton (determined by the 
FGS in September and adjusted for Triton’s different distance) was 
subtracted, and the remainder was divided by the flux from the star 
(also determined by the FGS in September) so that the full range of 
stellar flux corresponded to values between 0.0 and 1 .0. In the light- 
curve model fits, the zero level was fixed, but the full-scale signal 
from the star was a free parameter. The difference of the fitted values 
from 1 .0 shows that our calibration error is only a few tenths of one 
per cent. 

From fitting the entire light curve, the closest-approach distance 
between the centre of Triton’s shadow and the HST was determined 
to be 224 ± 4 km (first column of results in Table 2). This value 
places the shadow somewhat further north than predicted, but it is 
consistent with no detectable occultation at our Oahu station 
(Table 1). Using the closest-approach distance determined from 
the entire light curve, we fitted the main immersion and emersion 
sections of the light curve both separately and together (next three 
columns of results in Table 2). As a test of the self-consistency of our 
light-curve models, we fixed the atmospheric model parameters 
(“half-light radius”, “lambda at half-light,” and the “thermal- 
gradient exponent”; see ref. 8) at their values determined from 
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Triton, Neptune’s largest moon, has been predicted to undergo 
significant seasonal changes that would reveal themselves as 
changes in its mean frost temperature 1-3 . But whether this 
temperature should at the present time be increasing, decreasing 
or constant depends on a number of parameters (such as the 
thermal properties of the surface, and frost migration patterns) 
that are unknown. Here we report observations of a recent stellar 
occultation by Triton which, when combined with earlier results, 
show that Triton has undergone a period of global warming since 
1989. Our most conservative estimates of the rate of temperature 
and surface-pressure increase during this period imply that the 
atmosphere is doubling in bulk every 10 years — significantly 
faster than predicted by any published frost model for Triton 2 3 . 
Our result suggests that permanent polar caps on Triton play a 
dominant role in regulating seasonal atmospheric changes. Simi- 
lar processes should also be active on Pluto. 

The 4 November 1997 occultation of the star Trl80 (also known 
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Figure 1 Triton occultation light curves from the HST. Data before and after the 
occultation were used to establish the modulation of the signal due to the 
astrometer scan, which was then removed from the entire data set. The zero and 
full-scale stellar flux levels were established with photometric data from an earlier 
FGS visit to these objects on 11 September 1997. At that time, Triton’s magnitude 
as observed by the FGS was 13.4 and the magnitude of Tr180 was 10.6, The FGS 
data have been averaged at 1.0 s; a, immersion data (filled circles) and light-curve 
model fit (line); b, emersion data (filled circles) and light-curve model fit (line; see 
fits in Table 2). The zero point of the abscissa is arbitrary. The light-curve model 8 
used a power-law thermal gradient; residuals from the fit in a are shown in c, and 
the residuals from the fit in b are shown in d. The r.m.s. residual along the fuff 
signal is -0.0022 for a 1-s average, and is the result of photon noise. The 
remaining residuals that occur when the star is partially occulted are the result 
of unmodelled structure in Triton’s atmosphere and rarely exceed 0.01. The effect 
of these on our determination of the 1,400-km pressure can be estimated by 
differences among the fits in Table 1. The omitted central portion of the light curve 
(see text) corresponds to radii <1,400 km (altitudes <48 km); the atmosphere 
below that does not affect the 1,400-km pressure. 
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the combined immersion and emersion fit (second to last column in 
Table 2), and then fit the entire light curve again (last column in 
Table 2). 

Although the agreement is not perfect between all model fits for 
the atmospheric parameters at a radius of 1,400 km, it is very good. 
In particular, the close agreement between the immersion and 
emersion atmospheric parameters for two widely separated loca- 
tions on Triton supports the idea that the sublimation and con- 
densation of nitrogen maintains the surface frost at the same 
temperature. The pressure at 1,400 km, 2.3 ±0.1 |xbar, is derived 
from the fit to the immersion and emersion data (the error bar 
includes the systematic differences between the fits in Table 1). This 
is significantly greater than the 1,400-km pressure of 1.4 ± 0.1 |±bar 
measured with stellar occultations in 1995 9 and the value of 
0.8 ±0.1 juibar extrapolated from a surface pressure of 
14 ± t fxbar measured by Voyager in 1989 ,0,n . The pressure at 
1,400 km has certainly increased between 1995 and 1997, based 
on the two sets of stellar occultation data; the temperature at 
1,400 km has also increased from 47 ± 1 K (ref. 9) to 
50.3 ± 0.5 K. As these temperatures are consistent with the 1,400- 
km temperature predicted by atmospheric models based on Voyager 
data 12 ' 13 , a surface-pressure increase between the 1989 and more 
recent measurements is the most likely explanation for the differ- 
ence (rather than using an inappropriate model for extrapolating 
the pressure from the surface to 1,400 km). 


This surface-pressure increase implies a temperature increase of 
the surface frost, as the principal constituent of the tenuous atmos- 
phere, nitrogen, is presumed to be in vapour-pressure equilibrium 
with the surface frost. The surface pressures and corresponding 
equilibrium temperatures for N 2 are plotted in Fig. 2, where two 
cases are shown. 

In the first case, we have extrapolated the occultation pressures 
from a radius of 1,400 km to the surface radius of 1,352 km (an 
average of published values 14,15 ) using a multiplicative factor of 17.5, 
which is valid under the assumption that the shape^of the thermal 
profile 12 has not significantly changed since the Voyager encounter 
(filled symbols in Fig. 2). We justify this assumption by the 
consistency of the atmospheric temperature at 1,400 km (Table 2) 
with the models 12,13 . These measurements show a steady increase in 
pressure and surface-frost temperature. At this rate the atmosphere 
has tripled in bulk since the time of the Voyager encounter, which 
would require —0.9% of the total solar energy incident on Triton 
during the intervening time to sublime the frost. In the second case, 
we make a more conservative estimate of the surface pressure 
derived from the stellar occultation data by assuming an isothermal 
atmosphere below 1,400 km (open symbols in Fig. 2). Even for this 
case, the atmosphere has nearly doubled in bulk since the time of 
Voyager. 

Triton’s atmosphere is supported by surface frosts, which at any 
instant are expected to be at the same temperature everywhere, due 


Table 1 Observations 


Site 

Telescope aperture 
(m) 

Detector 

Filter 

Recording interval 
(uTon 4 Nov. 1997) 

Integration time 
(s) 

Aperture 

{arcsec) 

S/N* 

HSTt 

2.5 

4 PMTst 

F583WII 

04:59:32-05:24:44 

0.025 

5x5 

630 

LCC§ 

0.3 

CCD 

Nonell 

05:10:55-05:27:35 

0.125 

15 x 15 

33 


*S/N refers to" the occultation light curve {when the star is unocculted) over a time interval that corresponds to 20 km {approximately a scale height) of relative 

motion between the telescope and Triton's shadow. ... . . _ , 

t HST was located {in J2000 coordinates, kmrelative to the centre of the Earth) at (x.y.z) = (-3.906, -5,752, 761) for immersion, and (x T y,z) = (- 1 ,898, -6,484, 1,806) for emersion. Centre- 
line observations with NASA's Lear Jet Observatory were attempted, but failed due to target acquisition and tracking problems. Simultaneous optical and infrared observations with NASAs 
Infrared Telescope Facility on Mauna Kea (MKO) and optical observations with a 0.4-m portable telescope on Tern Island were foiled by clouds; other groups had similar bad weather on the 
University of Hawaii 2.2-m telescope on MKO and on Kauai {more details are provided at http://occult.mit.edu/trl80plans/). r 

t The four photomultiplier tubes (PMTs) of Astrometer 3 of the Fine Guidance Sensors were used in “TRANS” mode. We generated the light curve by adding the signals from all four PMTsn 
§The telescope at Leeward Community College (LCC) is located at longitude -157° 59' 6*E, latitude 2r23'30"N, and an altitude of 10 m. The data are a continuous series of CCD frames 
{each 60 by 57 pixels). The light curve was generated from these frames with synthetic aperture photometry of the blended image of Triton and the star. This station recorded data through the 
time of the event but saw no obvious occultation, which indicates that Triton's shadow path went somewhat north of the prediction (provided at http://occult.mit.eduAr180.html). 

HThe effective wavelength of the F583W filter was 583 nm; the open CCD used at LCC had an effective wavelength of 670 nm. 


Table 2 Models for the Triton occultation light curve obtained with the HST ____ 

Data selection 



All 

Immersion* 

Emersion* 

Immer. & emer. 


AI! 

Model parameters 
Background level 

0.00 

0.00 

0.00 

0.00 


0.00 

Background slope (10 -7 km -1 ) 

5.3±0.7 

0.00 

0.00 

2.2±0.8 


4.9±0.6 

Star signal 

0.9968 ±0.0003 

0.9957+0.0004 

0.9974 ±0.0003 

0.9966± 0.0002 


0.9971 ±0.0003 

Closest approach timet 

3.828 ±0.034 

3.808 

3.808 

3.808±0.015 


3.817+0.034 

Half-light radius (km) 

1,456.3+0.6 

1,456.3±0.3 

1,456.4+0.2 

1,456.3±0.2 


1,456.3 

Lambda (isothermal) at half-light 

65.3 ±0.4 

65.5±0.4 

66.4±0.3 

67.4 ±0.2 


67.4 

Thermal-gradient exponent 

1.9±0.3 

1.7±0.5 

4.4±0.4 

2.9±0.3 


2.9 

Minimum centre distance (km) 

224 ±4t 

224 

224 

224 


219-2 

Derived quantities atr = 1,400km (altitude 
Number density (10 u cm -3 ) 

= 48 km) 

3.23+0.09 

3.23±0.06 

3.42 ±0.04 

3.31 ±0.04 

] 

3.31 

Pressure (j^bar) 

2.25+0.06 

2.24±0.04 

2,36±0.03 

2.30±0.03 


2.30 

Temperature (K) 

50.6±0.3 

50.4 ±0.3 

50.1 ±0.2 

50.3±0.2 

S 

50.3 

Temperature gradient (K km' 1 ) 

0.07-0.01 

0.06 ±0.02 

0.16+0.01 

0.10+0.01 

= 

0.10 

Other derived quantities 
Minimum altitude probed (km)§ 

<25 

46 

46 

46 


<26 

Surface pressure (^bar) 

39±4 

39±4 

41+4 

40±4 

X 

40 

Fit information 
Degrees of freedom 

1,504 

321 

321 

644 

§ 

1,507 

x 2 per degree of freedomll 

17.8 

8.94 

4.22 

6.77 


18.1 


These models are for a pure N 2 atmosphere. Model parameters with error bars were fitted; those without were held constant. 

* The immersion data included the interval from 05:01 :03to 05:06:27 ut on 4 Nov. 1997; the emersion interval was 05:10:06 to 05: 15:30 ut. Immersion probed the atmosphere above -5.1 latitude 
and 10.4° longitude on Triton; emersion probed the atmosphere above -8.4° latitude and 206.4° longitude, 
t Seconds after 05:08:00 uron 4 Nov. 1997. p , , , . 

t The minimum centre distance is primarily controlled by the amplitude of the central flash, which will be discussed elsewhere. 

$ For the fits to all the data, the ‘minimum altitude probed’ refers to the light-curve centre. A surface radius of 1,352 km was used . 

\\ x 2 per degree of freedom is much greater than 1.0 for all fits, indicating that the deviations of the atmospheric structure from a power-law thermal gradient" {rather than photon noise) are the 
primary source of error {note the behaviour of the residuals in the lower panels of Fig. 1 ). 
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to efficient transfer of latent heat through the atmosphere 16,17 . This 
global frost temperature determines atmospheric pressure through 
vapour-pressure equilibrium; thus a frost temperature increase 
from 37.5 K (in 1989) to 39.3 K (in 1997) can be inferred from 
the observed increase in atmospheric pressure (the non-isothermal 
case in Fig. 2). 

The mechanisms that can cause a change in frost temperature 
(and hence surface pressure) are: (1) non-static (migrating) surface- 
frost distribution; (2) changes in the optical properties of the surface 
frost; and (3) changing insolation on a static surface-frost distribu- 
tion. (Although subsurface heat flow — geothermal or seasonal — 
can affect the energy balance, it is unlikely to have varied rapidly in 
the past decade.) 

Simple seasonal models of long-term frost migration, with uni- 
form frost albedo and emissivity 1,3 , show monotonic long-term 
transfer of most of Triton's nitrogen frost to the poles, which have 
the minimum seasonally averaged insolation. In most of these 
models, the remaining seasonal frost migrates between hemispheres 
in response to the seasonal insolation cycle, and all frosts should be 
almost exhausted from the southern hemisphere during the current 
extreme southern summer. This would result in a rapid decrease in 
atmospheric pressure, the opposite of what our 1995 and 1997 
observations indicate has occurred. 

A substantial decrease of mean frost albedo or emissivity could 
account for the pressure change. For instance, model "F 17 of Spencer 
and Moore 3 , which assumes unit emissivity, can account for the 
atmospheric change if the mean frost Bond albedo dropped by 0.12 
between 1989 and 1997. Observational 18,19 and theoretical 20 21 work 
suggests that large changes in frost optical properties may indeed 
occur on Triton. 

The frost’s total absorbed insolation may have increased, due to 
the increase in the southerly subsolar latitude in the past decade. We 
computed the change in surface pressure that would result solely 
from the change in subsolar latitude from August 1989 (47.8° south) 
to November 1997 (51.0° south), under the assumption that no 
volatile transport had occurred in the intervening 8 years. Five 
different maps of the 1989 N 2 -ice distribution on Triton were 
generated based on the albedos measured by Voyager. The emissiv- 
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Figure 2 Triton's global warming. The surface pressure and the corresponding 
vapour-pressure equilibrium temperature is plotted versus year of measurement. 
The measurement in 1989 refers to the results of an atmospheric model ,? based 
on several measurements made by Voyager 2. The 1995 data are from stellar 
occultations of Tr148A and Tr148B 9 , and the 1997 data refer to the present work. 
Each occultation datum has been plotted twice along with a corresponding least- 
squares fit: assuming that the atmosphere follows the model of Strobe! et al.' 2 
(filled circles, solid line), and making the conservative assumption that the 
atmosphere is isothermal below 1,400 km (open circles, dashed line). For the 
conservative assumption, the pressure increase is 1.4 ± 0.4 juLbaryr - ’; it is 
2.8 ± 0.8 p-baryr - ’ for the other case. In both cases, the pressure is seen to 
steadily increase, indicating that Triton has been undergoing a period of global 
warming. The distance of the points from the lines is larger than their error bars, 
which indicates that the warming of the surface frost may not have occurred 
linearly with time. 


ity of the N 2 ice was taken as a free parameter in each of these frost 
models, and was adjusted such that the temperature of the N 2 was 
equal to 37.5 K (corresponding to 14 jxbar surface pressure) in 1989. 
The N 2 -ice thermal balance was then recomputed using the subsolar 
latitude appropriate for Triton in late 1997, under the assumption 
that the distribution of N 2 ice was unchanged (albedos did change 
slightly, due to the illumination-angle dependence of the bolometric 
albedo, and the rotation of some N 2 - ice -covered portions of the 
globe into darkness). Several of these models result in pressure 
increases comparable to the measured increase. Thus, the observed 
global warming of Triton may be due to increased insolation of a 
permanent south polar cap. □ 
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